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HYPERFINE STRUCTURE OF THE GROUND 
STATE MUONIC 3 He ATOM 

A.A.Krutov, A.P.Martynenko 

Samara State University, Pavlov street 1, 443011, Samara, Russia 

On the basis of the perturbation theory in the fine structure constant a and the 
ratio of the electron to muon masses we calculate one-loop vacuum polarization and 
electron vertex corrections and the nuclear structure corrections to the hyperfine 
splitting of the ground state of muonic helium atom (pi e 2-ffe). We obtain total 
result for the ground state hyperfine splitting Av h f s = 4166.471 MHz which im- 
proves the previous calculation of Lakdawala and Mohr due to the account of new 
corrections of orders a 5 and a 6 . The remaining difference between our theoretical 
result and experimental value of the hyperfine splitting lies in the range of theoreti- 

i-Q ' cal and experimental errors and requires the subsequent investigation of higher order 

]> i 
i ■ corrections. 
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INTRODUCTION 



Muonic helium atom (// e %He) represents the simple three-body atomic system consisting 
of one electron, negative charged muon and positive charged helion i^He). Contrary to the 
muonic helium atom (/x e \He) it has more complicated ground state hyperfine structure 
which appears due to the interaction of the magnetic moments of the electron, muon and 
helion |l|-|9J]. The investigation of the energy spectrum of this three-particle bound state 
is important for the further check of quantum electrodynamics. Moreover, light muonic 
atoms (muonic hydrogen, muonic helium, ions of muonic helium etc.) represent a unique 
laboratory for precise determination of the nuclear properties such as the nuclear charge 
radius [10|. Hyperfine splitting (HFS) of the ground state of muonic helium atom (ji e \He) 



was measured many years ago with sufficiently high accuracy [11] : 



AzAf„ s = 4166.3(2) MHz. (1) 



There are two approaches to the calculation of the energy spectrum of muonic helium atom 
(/i e \He). First approach in |l|, |2|, \J\ is based on the perturbation theory (PT) for the 
Schrodinger equation. In this case there is the analytical solution for the three particle 
bound state wave function in the initial approximation. Using it the analytical calculation 
of different corrections to HFS can be performed. Contrary to the energy levels of two- 
particle bound states which were accurately calculated in quantum electrodynamics [12H17 



the hyperfine splitting of the ground state in muonic helium atom was calculated on the 
basis of the perturbation theory with essentially less accuracy. Another approach is built 
on the variatinal method [4j46|, |8|, |9|, |l8f|20| which allows to increase the accuracy of the 



calculation. In the beginning, the accuracy of the HFS calculation was not sufficiently high 
because corrections of six order in a were estimated only approximately. A feature that 



distinguishes light muonic atoms among the simplest atoms is that the structure of their 
energy levels strongly depends on the vacuum polarization, nuclear structure and recoil 
effects. Subsequently, the corrections of order a 2 to hyperfme splitting were studied in 
[5|, la] on the basis of variational and global-operator method. The theoretical value of HFS 
obtained in p, |6] contains very small uncertainty and agrees with the experimental value 

(1) - n 

In this work, which continues our investigation [21[, we aim to refine the calculation of 
Lakdawala and Mohr [l| using their analytical approach to the description of the muonic 
helium atom. We investigate such contributions of the one- loop electron vacuum polarization 
of order a^Me/M^ and the nuclear structure of orders a 5 , a 6 which are significant for the 
improvement of the theoretical value of the hyperfme splitting obtained in [l| on the basis of 
perturbation theory. Another purpose of our study consists in the improved calculation of 
the electron one-loop vertex corrections to HFS of order a 5 using the analytical expressions 
of the Dirac and Pauli form factors of the electron. 

The bound particles in muonic helium atom have different masses m e <C m^ <C m a . As 
a result the muon and helion compose the pseudo nucleus (fi \He) + and the muonic helium 
atom looks as a two-particle system in the first approximation. Three-particle bound system 
(/i e \He) is described by the Hamiltonian [22j, l23l ]: 



1 V 2 ^ V 2 _ 2a _ o_ 

2M U " 2M e e x, x e 



H = H + AH + AH rec , H = --—Vl--—V 2 e , (2) 



AH = — --, AH rec = -— V M • V e , (3) 

•^ lie *^e '"h 

where x M and x e are the coordinates of the muon and electron relative to the helium nucleus, 
M e = m e m h /{m e + m h ), M^ = m^m h /{m^ + m h ) are the reduced masses of subsystems 
[e \He) + and (a ^He) + . In the initial approximation the wave function of the ground state 



has the form j^-H]: 



*o(xe,x M ) = ^ e (x e )W(x M ) = -{2c?M e M,f l2 e- 2aM ^e- aM ^. (4) 

TV 

The hyperfme interaction in the ground state in (/x e \He) is determined by the following 
Hamiltonian: 

$H = -(fl N • /LI M )<J(X M ) - y (/V • Ve)5(x» - Xe) - y (>e • tlN)5(Xe), (5) 

where /x e = -g e e/(2m e )s e , ^ M = -g lx e/{2m l i)s^ Hn = -9Ne/{2m p )I N are magnetic mo- 
ments of the electron, muon and helion. The total spin of the three spin-1/2 particles can be 
either 3/2 and 1/2. The matrix element of (5) leads to the shift of the energy levels which 
takes on form: 

5E =< 5H >= -a (Ijv • s^) - b (s M • s e ) - c (s e ■ I N ), (6) 

where 

2ira g N g^ 2vra g^ 2vra g e g N 

a = — z- < d^j >, b = — ^ — < d(x M -Xe) >, c = — < d(Xe) > . (7) 

3 nipTrifj, 3 vn il m e 3 m e m p 



The diagonalization of the matrix element < SH > gives three eigenvalues: 

^1,2 = t(° + b + c) ± -(a 2 + b 2 + c 2 - ab - be - ca) 1/2 , z/ 3 = --(a + b + c). (8) 

The values z/ 12 and z/ 3 correspond to the total angular momentum | and |. In the case of 
muonic helium (/j, e \He) we have the relations a ^> b and a ^> c. So, the eigenvalues z/^2 
can be written with good accuracy in a more simple form: 

3 11 

v 1 = -a + ..., v 2 = --a + -(b + c) + .... (9) 

As a result the smaller hyperfme splitting interval related to the experiment (1) is given by 

Av h f° = v 2 -v z = l{b + c) + 0(~ ^). (10) 




Basic contributions to the coefficients a, b, c were calculated analytically in [1] from the 
contact interaction (5) in the first and second order PT. Taking into account numerical values 
of gyromagnetic factors g e = 2 for the b coefficient, g e = 2(l + /t e ) = 2(1 + 1. 1596521859- 10 -3 ) 
for the c coefficient, g^ = 2(1 + /^) = 2- (1 + 1.16592069(60) ■ 10^ 3 ), g N = 2-2.127497718(25), 
we obtain for them: 

8a(aM ) 3 

b = u F = — i ^- = 4516.307 MHz, b x = k„v f = 5.266 MHz, (11) 

Zrriem^ 

/;■.=--■/<,■( I-/.--,) i-3^ + -S 1/2 (^\ +(-f I ln-^ -( — ! S = -61. 322.1 IB:.. 

(12) 

Co = ^MKHhi = 1083.256 M#z, (13) 

4 m p 

Cl = Co (?*• + fMl) 2 ,,^ + (ln2 + I) (Ml) 2 ) = 8 .323 Mffz. (14) 

Note, that, as we determine contributions to the energy spectrum numerically, the cor- 
responding results are presented with an accuracy of 0.001 MHz. We express further 
the hyperfine splitting contributions in the frequency unit using the relation AE h ^ s = 
2irhAi> h f s . Modern numerical values of fundamental physical constants are taken from 
the paper [25]: the electron mass m e = 0.510998910(13) ■ 10~ 3 GeV, the muon mass 
m M = 0.1056583668(38) GeV, the fine structure constant a" 1 = 137.035999679(94), the 
helion mass ra(f-He) = 2.808391383(70) GeV, the muon anomalous magnetic moment k^ = 
1.16592069(60) • 10~ 3 , the electron anomalous magnetic moment n e = 1.15965218111(74) • 

io- 3 . 



II. EFFECTS OF THE VACUUM POLARIZATION 

The vacuum effects change the interaction (2), (3), (5) between particles in muonic helium 
atom. One of the most important contributions to HFS is determined by the one-loop 
vacuum polarization (VP) and electron vertex operator. Indeed, the vacuum loop leads to 





FIG. 1: The vacuum polarization effects. The dashed line represents the Coulomb photon. The 
wave line represents the hyperfine part of the Breit potential. G is the reduced Coulomb Green's 
function. 



additional factor a/ft in the interaction operator, so that corresponding correction to HFS 
is of the fifth order over fine structure constant. At the same time, the electron vacuum 
polarization and vertex corrections to the hyperfine splitting of the ground state contain the 
parameter equal to the ratio of the Compton wave length of the electron and the radius of 
the Bohr orbit in the subsystem ([i\He) + : m tl a/m e = 1.50886 .... It appears in the matrix 
elements with the use of the bound state wave function in which the characteristic momentum 
is of order m^a. It is impossible to use the expansion over a for such contributions to the 
energy spectrum. So, we calculate them performing the analytical or numerical integration 
over the particle coordinates and other parameters without an expansion in a. The effect 
of the electron vacuum polarization leads to the appearance of a number of additional 



corrections to the Coulomb potential which we present in the form [13|, |26 
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(16) 
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where x efl = |x e — x M |. They give contributions to the hyperfine splitting in the second 
order perturbation theory and are discussed below. In the first order perturbation theory 
the contribution of the vacuum polarization is connected with the modification of the hy- 
perfine splitting part of the Hamiltonian (5) (the diagram (a) in Fig.l). In the coordinate 



representation it is determined by the integral expression [27H29 
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(18) 
(19) 



Averaging the potential (18) over the wave function (4) we obtain the following contri- 
bution to the hyperfine splitting: 
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(20) 



There are two integrals over the muon and electron coordinates in Eq.(20) which can be 
calculated analytically: 

-^-^(x, _ Xe) = ; 2 - -3 , (2i) 
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They are divergent separately in the subsequent integration over the parameter £. But their 
sum is finite and can be written in the integral form: 
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(23) 

The order of this contribution is determined by two small parameters a and M e /M^ which 
are written explicitly. The correction hyp is of the fifth order in a and the first order in the 
ratio of the electron and muon masses. The contribution of the muon vacuum polarization to 
the hyperfine splitting is extremely small (~ 10 -6 MHz). One should expect that two-loop 
vacuum polarization contributions to the hyperfine structure are suppressed relative to the 
one-loop VP contribution by the factor a/n. This means that at present level of accuracy 
we can neglect these corrections because their numerical value is not exceeding 0.001 MHz. 
Higher orders of the perturbation theory which contain one-loop vacuum polarization and 



the Coulomb interaction (3) lead to the recoil corrections of order upotj^ In -^-. Such terms 
which can contribute 0.004 MHz are included in the theoretical error. 

Similar contribution to the coefficient c of order a 6 can be found analytically using the 
potential (19) («i = aM e /m e ): 



C yp = u F 



ag N m^ 



1 — «i(6ai + a\ — 3n) + (6 — 2>a\ + 6af) arccosai 



6nm r 



0.021 MHz. 
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Let us consider corrections of the electron vacuum polarization (15)-(17) in the second 
order perturbation theory (SOPT) (the diagram (b) in Fig.l). The contribution of the 
Coulomb electron-nucleus interaction (15) to the hyperfine splitting can be written as follows: 

167TOJ 
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(25) 



X E F T p 1 p e 2me ^ 1 ^0( X 2)V'e0(Xl), 

where the indices at the coefficient b indicate vacuum polarization contribution (VP) in the 
second order PT (SOPT) when the electron-helion Coulomb VP potential is considered. The 
summation in (25) is carried out over the complete system of the eigenstates of the electron 
and muon excluding the state with n, n! = 0. The computation of the expression (25) is 
simplified with the use of the orthogonality condition for the muon wave functions: 

32aM? f°° ,_. „ f°° o , f°° , - x , M^d+nisL) ~2x 3 (i+^- . 
bvp,soPT,e-h = v F ^ M t j i P(SWJ 4dx 3 J o x x dx x e xl M^+ aM J e { ™»J y 

(1 

(26) 
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= 0.150 MHz, 

where x < = min(xi,x 3 ), x > = max(xi,x 3 ), C = 0.577216... is the Euler's constant and 
Ei(x) is the exponential- integral function. It is necessary to emphasize that the transforma- 
tion of the expression (25) into (26) is carried out with the use of the compact representation 



for the electron reduced Coulomb Green's function obtained in Refs. [30 
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The contribution (26) has the same order of the magnitude 0(a 5 ^) as the previous cor- 
rection (23) in the first order perturbation theory. Similar calculation can be performed 
in the case of muon- nucleus Coulomb vacuum polarization potential (16). The intermedi- 
ate electron state is the IS state and the reduced Coulomb Green's function of the system 
appearing in the second order PT transforms to the Green's function of the muon. The 
correction of the operator (16) to the hyperfine splitting (coefficient b) is obtained in the 
following integral form: 



l | rn e £ 
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0.048 MHz. 



The vacuum polarization correction to HFS which is determined by the operator (17) in the 
second order perturbation theory is the most difficult for the calculation. Indeed, in this 
case we have to consider the intermediate excited states both for the muon and electron. 
Following Ref . [22J we have divided total contribution into two parts. The first part in which 
the intermediate muon is in the IS state can be written as: 



256a 2 {aM e ) 3 {2aM fl f >-x 
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a;i(16a 2 M 2 -4m 2 £ 2 ) 2 
After the substitution (30) in (29) the numerical integration gives the result 

b vp, soPT{n = 0) = -0.029 MHz. 



(31) 



Second part of the vacuum polarization correction to the hyperfme splitting due to the 
electron- muon interaction (17) can be presented as follows: 



b VP, SOPT, tte(n ^ 0) = 
X ^^n(x 3 )^„(x 2 ) 
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In the expression (32) we have replaced the exact electron Coulomb Green's function by 
the free electron Green's function which contains B = \2M P {E jm — E^ — E eG ] 1 / 2 . (see more 
detailed discussion of this approximation in Re fs.|22l .l23l|). We also replace the electron wave 
functions by their values at the origin as in Ref. [22J neglecting higher order recoil corrections. 
After that the integration over xi can be done analytically: 
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where we have performed the expansion of the first exponential in the square brackets over 
powers of £>|x 3 — x 2 |. As discussed in Ref. |22| one can treat this series as an expansion over 



the recoil parameter J M e /M^. For the further transformation the completeness condition 
is useful: 

2 ^n(x 3 )^* n (x 2 ) = <5(x 3 - x 2 ) - ^o(x 3 )^ (x 2 ). (34) 

The wave function orthogonality leads to the zero results for the second and fifth terms in 
the square brackets of Eq.(33). The first term in Eq.(33) gives the leading order contribution 
in two small parameters a and M e /M fJb : 
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3a 2 M P 



8m P 



-u F , 



(35) 



'12 



V F 



2a 2 

6n M e 



p(0 



d£ 



M?,a A 



£ (AaM^ + 2m e £Y 



256 + 232^ + 80^ + 10 <? 



M^a 



Mia 2 



M 3 a 3 



(36) 



The numerical value of the sum 611+612 = —0.062 MHz is included in Table I. It is important 
to calculate also the contributions of other terms of the expression (33) to the hyperfine 
splitting. Taking the fourth term in Eq.(33) which is proportional to B 2 = 2M e (E fJ-n — E^q) 
we have performed the sequence of the transformations in the coordinate representation: 
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Evidently, we have the divergent expression in Eq.(37) due to the presence of the ^-function. 
The same divergence occurs in the other term containing B 2 in the square brackets of Eq.(33). 
But their sum is finite and can be calculated analytically with the result: 
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Numerical value of this correction 0.009 MHz is smaller than the leading order term. Let 
us consider also the nonzero term in Eq.(33) proportional to B 3 . First of all, it can be 
transformed to the following expression after the integration over £: 
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where the sum S 3 / 2 is defined as follows: 



s P = E 



E 



fin 



E 



fiO 



R,, 



< 4>fio\ — \lpim 

a. 



> 



(40) 



R„ = 2a 2 M, 



fn 



\f2aM il . Using the known analytical expressions for the dipole matrix 



r~ r" ' r" ' r" v -' ** ■*■ *■ 1 ■ 1 

elements entering in Eq.(40) in the case of the discrete and continuous spectrum 12|, [31 
can write separately their contributions to the sum S3/2 in the form 
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As a result S3/2 = 3.2722 .... The similar calculation of the sum S1/2 relating to this 
problem (see Ref.|22[) gives Si/ 2 = 2.9380 . . .. Numerical value (39) is taken into account in 
the total result presented in Table I. 

The Coulomb vacuum polarization (15) does not contain the muon coordinate, so, its 
contribution to the coefficient c in the second order PT can be derived taking n = for the 




FIG. 2: Vacuum polarization effects in the second order perturbation theory. The dashed line 
represents the first part of the potential AH (3). The wave line represents the hyperfine part of 
the Breit potential. 



muon state in the Coulomb Green's function. Moreover, the <5(x e ) function in (5) results in 
the appearance of the electron Green's function with one zero argument: 
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-0.044 MHz. 



The vacuum polarization in the Coulomb {jx — N) interaction does not contribute to c in 
SOPT because of the orthogonality of the muon wave functions. Let consider correction to 
the coefficient c arose from (17) in SOPT. Only intermediate muon state with n = in the 
Green's function gives the contribution in this case. Using (43) we make integration over 
electron coordinates and present this correction in the form (7 = m e £,/2aM fJj , 7 X = M e /4M At ): 
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[9+27 1 (37+47 1 (17+97i))]+47 1 (l+27 1 )(7+27 1 )(l+7 2 +27(l+27 1 )+27 1 (3+27 1 ))ln(2 7l )] = 

= 0.009 MHz. 

There exist another contributions of the second order perturbation theory in which we 
have the vacuum polarization perturbation connected with the hyperfine splitting parts 
of the Breit potential (18), (19) (see Fig. 2). Other perturbation potential in this case is 
determined by the first term of relation (3). We can divide the HFS correction of (18) into 
two parts. One part with n = corresponds to the ground state muon. The other part with 



10 

n 7^ accounts the excited muon states. The ^-function term in Eq.(18) gives the following 
contribution to HFS at n = (compare with Ref.|22|): 



a r°° HM e 
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Obviously, this integral in the variable £ is divergent. So, we have to consider the contribution 
of the second term of the potential (18) to the hyperfme splitting which is determined by 
the more complicated expression: 
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N , ^ T/w(x 3 )^* n ,( Xl ) 

X Z, p P AV 2 (x 1 )^ e0 (xi) 



-EeO — E en i 



where 



/g— 2m e f|x3— X4| 
dx 4 C (x 4 )- j-^/io(x4) = (48) 

Xa — X4 



4(2aM p 



\3 



, ( , .SaM.e"^^ + e" 4 "^ (-8aM, - 16a 2 M 2 x 3 + 4m 2 £ 2 x 3 ) 

x 3 [(4aM M ) i - (2m e £) 2 J 2 L p v p ^ ' 

Ay 2 ( Xl ) = /dx 2 Wo(x 2 ) f , " 1 - -) ^o(x 2 ) = --(1 + 2aM^ 1 )e- 4aM ^ 1 . (49) 

J \|X 2 -Xi| Xi/ Xi 

Nevertheless, integrating in Eq.(47) over all coordinates we obtain the following result in 
the leading order in the ratio {M e /M IJ ): 

, , m m e Ml r ^^^ 32 + 63 7 + 44 7 2 + ll7 3 

6 ^, sew, 12 (n = 0) = -^j^gg^ / p(£« ^^ • (50) 

This integral also has the divergence at large values of the parameter £. But the sum of 
integrals (46) and (50) is finite: 

aM P r°° 11 + 12t + 3^ 2 

b vp, soft, 11 (n = 0)+b VPt soft, 12 n = 0) = vf-^-t \ p(£)d£ .,',, ! = 0.008 M#z. 

487rM M A (I + 7) 4 

(51) 

Let us consider now the terms in the coefficient b with n^0. The delta-like term of the 

potential (18) gives the contribution to the HFS known from the calculation of Ref.[22]: 

/ , „s a f°° ,_ x ,_ / 35M e \ ,„„, 

»k,.»,™( B *<>) = ,,— / ^--^j. (52) 

Another correction from the second term of the expression (18) can be simplified after the 
replacement the exact electron Green's function by the free electron Green's function: 



b vp, soft, 22{n ± 0) = - e —± / p{i)edi / rfx 2 / rfx 3 x (53) 

/ e -2m e £|x3-x 4 | 00 
dx 4 V>*o(x 4 ) — — ^ ^( X 4)^( X 2)I X 3 - X 2 |-0 M o(x 2 ) 
Xa — Xi /n 



The analytical integration in Eq.(53) over all coordinates leads to the result: 
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b vp, soft, 22 {n 7^ 0) — ~~ U F 



aM e 
3vrM„ 



p(m 



;i+7)^ 



f . 1 in 215 7 2 35 7 4 



(54) 
The sum of expressions (52) and (54) gives again the finite contribution to the hyperfine 
splitting: 

b vp, sopt, 21 (n 7^ 0) + Az/J£ 50PT 22 (n 7^ 0) = (55) 



-u F 



aM e 
3ttM ; 



p(0# 



35 + 767 + 597 2 + I67 3 



-0.062 MHz. 



16(1 + 7 ) 4 

Despite the fact that the absolute values of the calculated VP corrections (31), (35), (36), 
(37), (39), (51), (55) are sufficiently large , their summary contribution to the hyperfine 
splitting (see Table I) is small because they have different signs. 

The hyperfine splitting interaction (19) gives the contributions to the coefficient c in 
second order PT. Since the muon coordinate does not enter into the expression (19), we 
should set n = for the muon intermediate states in the Green's function. The initial 
formula for this correction is 



c VP, SOPT 



8a 3 gN 
9TTm P m r 



P{£)d£ j dXl/ dx 3 y rfx 4 |^o (X 3 ) | 2 V'eo( X 4)V'eo(Xl ) X (56) 



X 



x 3 -x 4 



1 

X i 



G e (x 4 ,Xi) 7t5(x! 






The integration over x 3 can be done analytically as in (49). Then it is useful to divide 
(56) into two parts. The coordinate integration in the first term with the 8 - function is 
performed by means of (43). In the second term of (56) we use the electron Green's function 
in the form (27). The summary result can be presented in the integral form in the leading 
order in M e /M^. 



C VP, SOPT 



2a 5 g N M ( 



97rm e m p M At 



p(£K 



°^ Z 2aM u 



1 + 



2aM u 



0.013 MHz. 



(57) 



III. NUCLEAR STRUCTURE AND RECOIL EFFECTS 



Another set of significant corrections to the hyperfine splitting of muonic helium atom 
which we study in this work is determined by the nuclear structure and recoil [l4j, 132^34 . 
The charge and magnetic moment distributions of the helion are described by two form 
factors GE{k 2 ) and Gm(^ 2 ) for which we use the known parameterization [35j: 



G E (k 2 



b 2 k 2 e 



+ d 



(fc+jo)_ 



(fc-jpr 



(58) 



G M (k 2 



e"^ 2 



b\k 2 e 



(59) 
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where numerous parameters a, b, c, 5i, &i, ci, go, p are written explicitly in [35]. In I7 
interaction the nuclear structure correction to the coefficient c is determined by the ampli 
tudes shown in Fig. 3. Purely point contribution in Fig.3(b) leads to the HFS value (13) 
Then the nuclear structure correction is given by 

(1 + K^gtfmn 



C str, I7 — Vp- 



2m rj 



[G M {x)e- 2aM ° x dx-l 



-0.072 MHz. 



(60) 



FIG. 3: Nuclear structure correction to coefficient c in I7 interaction. The bold point represents 
the nuclear vertex operator. The wave line represents the hyperfine part of the Breit potential. 



+ 




G-f 



FIG. 4: Nuclear structure corrections to coefficient c in 27 interactions. The bold point represents 
the nuclear vertex operator. The wave line represents the hyperfine part of the Breit potential. 
Dash line corresponds to the Coulomb potential. 



G 



-6- 




G 



G 



FIG. 5: Nuclear structure corrections to coefficient c in second order PT. The bold point represents 
the nuclear vertex operator. The wave line represents the hyperfine part of the Breit potential. 
Dash line corresponds to the Coulomb potential. G is the reduced Coulomb Green's function. 

Two-photon (e — h) interaction amplitudes (see Fig. 4) give the contribution to HFS of 
order a 5 . It can be presented in terms of the Dirac and Pauli form factors F± and F2 13, [28 



J 2 « 5 



'Ml 



C str, 27 



3iTm e mh 



no 



fv-rn, 



(61) 



13 



V(k) 



2F 2 k 2 



+ 



M P 



rn, 



Wh (m e - m h )k(k + J4m 2 e + k 2 ) 



„„,» ^o,, 32F 2 m 2 k 2 AF 2 k 4 AF 2 k A 
+64F 1 F 2 k 2 + lQF£k 2 + 2 e + ' 



-128^77^ - 128FiF 2 m^ + 16F 2 k 2 + 
M P . 



mi 



mi 



mi 



+ 



(m e — mh)k(k + y 4m| + k 2 ) 



x 



x 



\2%Flm 2 h + 128FiF 2 m^ - lQF 2 k 2 - ^F x F 2 k 2 - 48F 2 k 2 
The subtraction term in Fig.4 is taken as follows: 

64M*Z 2 a 5 g N r°° dk 



C iter, str 



:G M (k 2 



(62) 



3 6m e m/ l 7r Jo k 2 

It is only the part of the iteration contribution < Vi 7 xG'x V\ 1 > s {r ■ Other parts are used 
in the second order PT (see Fig. 5). As a result we obtain: 

C s tr, 2 7 + C iter, str = "0.077 MHz. (63) 



G 



FIG. 6: Nuclear structure correction to coefficient b in the second order perturbation theory. The 
wave line represents the hyperfine (e — fj,) interaction. G is the reduced Coulomb Green's function. 

The nuclear structure corrections to the coefficient c in second order PT which are pre- 
sented in Fig. 5. Here we have two different contributions. First contribution is related with 
diagrams in Fig. 5 (a), (b) when the hyperfine part of one perturbation is determined by mag- 
netic form factor Gm and the other perturbation is described by the charge radius of the 
helion r^. 

AV£__ h (r) = 1-nZa < r\ > 6(r). (64) 



The general integral structure of this correction and its numerical value are the following: 

_ h = Vp A7{Z 9Nm^ < „ r xdxe -2xt 4x n n2x + C ^ +4x 2_ 1Q , Gi 

QaM e mv Jo 



Cl, str, SOPT,e- 



L e" L p 



aM P 



(65) 



0.00005 MHz. 



Numerical value of the contribution c\ t s t r _ sopr,e-h is obtained by means of the charge 
radius of the helion r^ = 1.9642(11) fm 36]. The second nuclear structure contribution 
from the diagrams Fig. 5 (c), (d) is evaluated by means of the potential AH (3) and the 
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helion magnetic form factor. For the amplitude in Fig. 5(c) we make the integration over the 
muon coordinate in the muon state at n = and present the correction to the coefficient c 
in the form: 

c 2 , sir, soPT,e-h+ci = — — — — / x e 2M " G M {——)dx / y(l+-)e 2M »'dyx 

24m e m p M^ Jo AaM^ Jo 2 



(66) 
M P s , , M e „., M e , 7 _ M P , i- e 7i% x < 



M, 



ln(-^-x>) - ln(— ^-x<) + Ei(— ^x<) + - - 2C - -^-{x + y) + 



^*> ■ V 2M, >; V 2M, <; V 2M M <; 2 2M/ y; ^-x< 

= 8.249 MHz. 

Subtracting the point contribution c\ (14) we find c<i. s tr, sopr,e-h — —0.074 MHz. 

There is the nuclear structure contribution to the coefficient b in second order PT which 
is presented in Fig.6. If we consider the Coulomb interaction between the muon and helion, 
then the structure correction takes on form: 



64vr 2 a 2 

'str, ll— h 7> 

9m e m M 



r 2 ^=(2aM / ,) 3/2 |dx3^ (x3)|Veo(x 3 )rG'^X3,0,^ ). (67) 



After that the analytical integration over the coordinate x 3 in Eq.(67) can be carried out 
using the representation of the muon Green's function similar to expression (43). The result 
of the integration of order 0(a 6 ) is written as an expansion in the ratio M e jM^\ 

W-* = -v F \a 2 My a U^ - y ^f + ..) = -0.010 MHz. (68) 

The same approach can be used in the calculation of the electron-nucleus interaction. The 
electron feels as well the distribution of the helion electric charge. The corresponding con- 
tribution of the nuclear structure effect to the hyperfine splitting is determined by the 
expression: 

QAtt or r r 
b s tr,e-h = r\ J dxi / dx 3 |^* (x 3 )| 2 ^ e o(x3)G' At (x3,x 1 , J E e0 )^ e o(x 1 )(5(x 1 ). (69) 

Performing the analytical integration in Eq.(69) we obtain the following series: 
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J s tr,e-h = -v F ^-a 2 Mlr 2 a 



5 _ ln ^ + 4( 3 ln^-7U4(H-31n^ 
M, Ml \ M, ) Ml \ 2 M„, 



(70) 
= -0.003 MHz. 

We have included in Table I the total nuclear structure contribution to the coefficient b 
which is equal to the sum of the numerical values (68) and (70). 

Special attention has to be given to the recoil corrections connected with two-photon 
exchange diagrams shown in Fig. 7 in the case of the electron-muon interaction. The lead- 
ing order recoil contribution to the interaction operator between the muon and electron is 
determined by the expression |5j, [l3|, l37| : 



a 2 , m 



AKS-e(^) = -8 2 2 In ^(s M s e )<5(x Me ). (71) 



m 2 — ml a/ 
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a b 

FIG. 7: Two photon exchange amplitudes in the electron-muon hyperfine interaction. 

Averaging the potential AV rec s 4l _ e over the wave functions (4) we obtain the recoil correction 
to the coefficient b: 

Orec,a-e = "f 5 ^ In — ^ = 0.812 MHz. (72) 

it m — m z e m e 

There exist also the two-photon interactions between the bound particles of muonic helium 
atom when one hyperfine photon transfers the interaction from the electron to muon and 
another Coulomb photon from the electron to the nucleus (or from the muon to the nucleus). 
Supposing that these amplitudes give smaller contribution to the hyperfine splitting we 
included them in the theoretical error. 

IV. ELECTRON VERTEX CORRECTIONS 

In the initial approximation the potential of the hyperfine splitting is determined by 
Eq.(5). It leads to the energy splitting of order a 4 . In QED perturbation theory there is the 
electron vertex correction to the potential (5) which is defined by the diagram in Fig.8(a). 
In momentum representation the corresponding operator of hyperfine interaction has the 

form: 

"a 2 /cr cr ' 



where G^(k 2 ) is the electron magnetic form factor. We extracted for the convenience the 
factor a/ir from G^{k 2 ) — 1 . Usually used approximation for the electron magnetic form 

factor Gjtf(k 2 ) ~ G^i (0) = 1 + n e is not quite correct in this task. Indeed, characteristic 
momentum of the exchanged photon is k ~ aM^. It is impossible to neglect it in the 
magnetic form factor as compared with the electron mass m e . So, we should use exact 



one-loop expression for the magnetic form factor which was obtained by many authors [26 
Let us note that the Dirac form factor of the electron is dependent on the parameter of the 
infrared cutoff A. We take it in the form A = m e a using the prescription m e a 2 C A C m e 
from Ref . 12J . 



Using the Fourier transform of the potential (73) and averaging the obtained expression 
over wave functions (4) we represent the electron vertex correction to the hyperfine splitting 
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G 



G 



a b c 

FIG. 8: The electron vertex corrections. The dashed line represents the Coulomb photon. The 
wave line represents the hyperfine part of the Breit potential. G is the reduced Coulomb Green's 
function. 

as follows: 



'vertex, I7 




k 2 dk 






r M 



(*' 



(74) 



1 + 



m e 



4aM„ 



' M e 

2M„ 



+ 



m e 



4aM„ 



2} 



4.218 MHz. 



Let us remark that the contribution (74) is of order a 5 . Numerical value (74) is obtained 
after numerical integration with the one-loop expression of the electron magnetic form factor 



G^,{k 2 ). If we use the value G^(k 2 = 0) then the electron vertex correction is equal 
5.266 MHz. So, using the exact expression of the electron form factors in the one-loop 
approximation we observe the 1 MHz decrease of the vertex correction to the hyperfine 
splitting from I7 interaction. Taking the expression (73) as an additional perturbation 
potential we have to calculate its contribution to HFS in the second order perturbation 
theory (see the diagram in Fig. 8(b)). In this case the dashed line represents the Coulomb 
Hamiltonian AH (3). Following the method of the calculation formulated in previous section 
(see also Refs.[22[|23|) we divide again total contribution from the amplitude in Fig. 8(b) into 
two parts which correspond to the muon ground state (n = 0) and muon excited intermediate 
states (n ^ 0). In this way the first contribution with n = takes the form: 



b vertex, SOPT{n — 0) 



8a< 



3ir 2 m P m, 



k 



g m(^ 2 ) - 1 dk I dxi / dx 3 ^ e0 (x 3 ) 



x 



(75) 



x AV^k, x 3 )G e (xi, x 3 )AV 2 (x 1 )^ e0 (xi), 
where AV^Xi) is defined by Eq.(49) and 

at>/l n fj i / ^sin(fc|x 3 -x 4 |) sin (i^r) 

AVi(fc,X3) = / dx 4 Vvo(x 4 ) i i t/Vo(x 4 ) = 



|x 3 -x 4 | 



•i'3 



1 + 



fc 2 
(4aM M ) 2 



2 ' 



(76) 
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Substituting the electron Green's function (27) in Eq.(75) we transform desired relation to 
the integral form: 



b vertex, SOPT(n = 0) = Up 



a 



m P 



2 ( M P \ 2 rook 



16tt 2 \ aM u 



M„, 



( e )(U2\ 



Gtlik 



dk 



1 + 



I'yk 2 



(4aM,) 2 



(77) 



X / X-xt 

/o 



-x 3 



( m z k \ j f°° 

s ni^ 3 rWo xi 



(i+l) e -( i+ ^)& i: 



2M M 
M e x > 



ln( - ; ', u- 



, , M e , ,_, , M e . 7 ^ M e . . 1 - e 



2M U - L < 



2Mu< 



-0.208 M#z. 



One integration over the coordinate x\ is carried out analytically and two other integrations 
are performed numerically. Second part of the vertex contribution (Fig. 8(b)) with n ^ can 
be reduced to the following form after several simplifications which are discussed in section 
II (see also Refs.0 0): 



8a 4 M P M, 3 



^vertex, SOPT{n ^ 0) — Up — 

ir 6 



j e -2aM,x 2d ^ 2 f e -«M e x 3rfx3 f e -2«Af M x 4dx4 x ^ g j 



x / fcsin(fc|x 3 - x 4 |) (G^)(k 2 ) - l) — ^ [5(x 4 - x 2 ) - V>( x 4)VVo( x 2)] • 

/O v 7 IX3 — X4I 

We divide expression (78) into two parts as provided by two terms in the square brackets of 
(78). After that the integration (78) over the coordinates x l5 x 3 is carried out analytically. 
In the issue we obtain (72 = m e k/AaM^): 



W,vertex, SOPT^n ^ 0) — Vp 



ft 



m R 



32tt 2 \ftM„ 



M r°° 



M^jo 



47i(7? " 1) 7i(3 + 7?) , 4 7l 2 (7! 2 -l) , 



G { Sik 2 ) 
1 + H 



dk 



1 



1 + 12 



2\3 



(1 + 7S) 



2\2 



^2,-uertex, SOPT 



n/0) 



ft 



(7i 2 + 7l) 3 (7i 2 + 7 2 2 ) 2 
3 



(7i 2 " I) 3 
2.528 MHz 



x 



(79) 



m P 



-v F 



32tt 2 VftM, 



Me 

M M A, 



k 2 



1 + 7 2 2 



2N2 



(71 + 1) 



(7i 2 + 7l) [(l + 7i) 2 + 7 2 2 ] 2 (7i + l) 2 + 7l (7i 2 +7l) 



Gi?(A; 2 )-l 

7 2 2 ~ 3 7l 2 
2^3 



dk x 



(80) 



-0.831 MHz. 



It is necessary to emphasize that the theoretical error in the contributions b 12 , vertex, so <pr( n ^ 



0) is determined by the factor JM e jMn connected with the omitted terms of the expansion 
similar to Eq.(33) (see also Refs. |22l. |23_|). It can amount to 10% of the results (79), (80) 
that is the value near 0.2 MHz. 

Until now we consider the electron vertex corrections connected with the hyperfine part 
of the interaction Hamiltonian (5). But in the second order perturbation theory we should 
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analyze vertex corrections to the Coulomb interactions of the electron and muon, electron 
and nucleus. Then in the coordinate representation we have the following potential: 



rC 



£** vertex, eN\ x e) + £\ * vertex, eiApe-P) 



2a' 



7T- 



t(e) 



GW(k 2 



k 



■dk 



sin(fcx. 



eft / 



sin(kx e 



x 



e/i 



X r 



(81) 

where we extract again the factor a /it from G% (k 2 ) — 1 . G F is the electron electric 
form factor. One part of the contribution in Fig. 8(c) is specified by the electron- muon 
intermediate states in which the muon is in the ground state n — 0. This correction is 
determined by both terms in large parentheses of Eq.(81) and can be presented as follows: 



x / x\e 

/o 



bc,vertex, SOPr{n — 0) — 

\G { E\k 2 ) 



a 



v F - 



7T Z 



Me 
M„ 



XqG 



-»(i+^r) 



dxs x 



(82) 



2M M 
M e x > 



ln( 



M e 
2M, 



-xi 



■X, 



dx\ 



-ln( 



dk 



k 



sin 




lik' 2 



(4aAf M ) s 



+ 1 



M e 
2M, 



M P 



-a;- 



+ K <^<»4- 2C - ,u.. 



-1.303 MHz. 



.xi +a; 3 ) + 



1-e 



2M U , X < 



2M, 



x- 



The index " C" means that the vertex correction to the Coulomb part of the Hamiltonian 
is considered. Excited states of the muon (n 7^ 0) contribute to another part of the matrix 
element (Fig.8(c)). Changing the Coulomb Green's function of the electron by free Green's 
function (see discussion in section II) we can carry out the coordinate integration and express 
the correction to HFS as one-dimensional integral: 



JC,vertex, SOPT 



>i/0) 



8a M P faM, 



-v F 



vr 2 M„ 



m e . 



v(e) 



G^>(k 2 



dk 



1 



k 2 



i'ik z 



(4aM M ) 2 



(83) 



1.806 MHz. 



The electron vertex corrections investigated in this section have the order a 5 in the hyperfme 
interval. Summary value of all obtained contributions (74), (77), (79), (80), (82), (83) is 
equal to 6.210 MHz. It differs by a significant value 0.944 MHz from the result 5.266 MHz 
which was used previously by many authors for the estimation of the electron anomalous 
magnetic moment contribution. On our opinion, it is necessary to use the same approach 
for the calculation of the electron vertex corrections by the variational method 19( in which 
the bound state wave function has the form ^(x e ,x M , x e/i ) = J2i de~ aiX& ~^ iX ^~ liXe ^ and the 
parameters «;, Pi, 7, are chosen randomly between some minimal and maximal values. 



V. CONCLUSIONS 



In the present study, we have performed the analytical and numerical calculation of several 
important contributions to the hyperfme splitting of the ground state in muonic helium atom 
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TABLE I: Hyperfine splitting of the ground state in the muonic helium atom (// e zHe). 



Contribution to the HFS 


b, MHz 


c, MHz 


Reference 


The Fermi splitting 


4516.307 


1083.256 


(11), (13), [1,22,23] 


Recoil correction of order 

a 4 (m e /m M ) 


-64.322 


8.323 


(12), (14), [1] 


Correction of muon anomalous 
magnetic moment of order a 5 


5.266 





(11), [1,22] 


Recoil correction of order 


0.105 




[22,23] 


a 4 (M e /m Q ) A /(M e /M M ) 


Relativistic correction of order a 6 


0.040 


0.087 


[2,5] 


One-loop VP contribution in I7 
interaction of orders a 5 , a 6 


0.036 


0.021 


(23), (24) 


One-loop VP contribution in the 
second order PT 


0.062 


-0.023 


(26),(28),(31),(35),(36),(38), 
(39),(44),(45),(51),(55),(57) 


Nuclear structure correction 
in I7 interaction of order a 6 





-0.072 


(60) 


Nuclear structure and recoil correction 
in 27 interactions of order a 5 





-0.077 


(63) 


Nuclear structure correction of order 
a 6 in second order PT 


-0.013 


-0.074 


(14),(65),(66) 


Recoil correction of order 
a 5 (m e /m^) ln(m e /m M ) 


0.812 


— 


(72),[5,37J 


Electron vertex correction of order a 6 


-0.615 


-0.035 


[13, 38-40] 


Electron vertex contribution 
of order a 5 


6.210 


— 


(74),(77),(79), 
(80), (82) ,(83) 


Summary contribution 


4463.888 


1091.406 


Al/ hfs _ 4166 .4 71 MHz 



connected with the vacuum polarization, the nuclear structure, recoil effects and the electron 
vertex corrections. To solve this task we use the method of the perturbation theory which 
was formulated previously for the description of the muonic helium hyperfine splitting in 
Refs.[ll, |22|, |23J. We have considered corrections of the electron vacuum polarization, electron 
electromagnetic form factors and the nuclear structure effects of orders a 5 and a 6 . The 
numerical values of the corresponding contributions are displayed in Table I. We present 
in Table I the references to the calculations of other corrections which are not considered 
here. The relativistic correction was obtained in Ref . [3|, [5| , the electron vertex corrections of 
order a(Za)Ep were calculated in the case of hydrogenic atoms in Refs. 13|, I38H40J . Basic 
contributions to the hyperfine splitting obtained by Lakdawala and Mohr are also included 
in Table I because our calculation is closely related to their approach. 

Let us list basic points related to the calculation. 

1. For muonic helium atom, the vacuum polarization effects are important for obtaining 
the high accuracy of the calculation. They give rise to the modification of the two-particle 
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interaction potential which provides the a^-M^-order corrections to the hyperfine structure. 
The next to leading order vacuum polarization corrections (two-loop vacuum polarization) 
are negligible. 

2. The electron vertex corrections to the coefficient b should be considered with the 
exact account of the one-loop electromagnetic form factors of the electron because the char- 
acteristic momentum incoming in the electron vertex operator is of order of the electron 
mass. 

3. The nuclear structure corrections to the ground state hyperfine splitting are expressed 
in terms of electromagnetic form factors and the charge radius of the helion. 

4. Analyzing the one-loop electron vacuum polarization and vertex effects and the nu- 
clear structure contributions in each order in a, we have taken into account recoil terms 
proportional to the ratio of the electron and muon masses. 

The resulting numerical value 4466.471 MHz of the smaller ground state hyperfine split- 
ting in muonic helium (// e \He) is presented in Table I. It is sufficiently close both to the 
experimental result (1) and the earlier performed calculations by the perturbation theory, 
variational approach in |l|, |2|, |5j . The estimation of the theoretical uncertainty can be done 
in terms of the Fermi energy vp and small parameters a and the ratio of the particle masses. 
On our opinion, there exist several main sources of the theoretical errors. First of all, as we 
mentioned above comprehensive analytical and numerical calculation of recoil corrections of 
orders a^ffs a 4 T7t, °^-j^k ^(M^/M e ) was carried out by Lakdawala and Mohr in the sec- 



ond order PT in Refs. 22|, |23j. The error of their calculation connected with the correction 



u Fjpi l n jf- consists 0.6 MHz. The second source of the error is related to contributions of 

order o?vp fa 0.2 MHz which appear both from QED amplitudes and in higher orders of the 
perturbation theory. Another part of the theoretical error is determined by the two-photon 
three-body exchange amplitudes mentioned above. They are of the fifth order over a and 
contain the recoil parameter (m e /m a ) ln(m e /m a ), so that their possible numerical value can 
be equal ±0.05 MHz. Finally, a part of theoretical error is connected with our calculation 
of the electron vertex corrections of order a 5 in section IV. It consists at least 0.2 MHz 
(see the discussion after Eq.(80)). We neglect also the electron vertex contributions of order 
upaMe/M^ fa 0.2 MHz which appear in higher orders of the perturbation theory. Thereby, 
the total theoretical uncertainty is not exceeded ±0.7 MHz. The existing difference between 
the obtained theoretical result and experimental value of the hyperfine splitting (1) equal 
to 0.171 MHz lies in the range of total error. Theoretical error which remains sufficiently 
large in the comparison with the experimental uncertainty, initiates further theoretical in- 
vestigation of the higher order recoil contributions and more careful construction of the 
three-particle interaction operator connected with the multiphoton exchanges. 
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